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Abstract. The accuracy of Monte Carlo Glauber model descriptions of minimum- 
bias multiplicity frequency distributions is evaluated using data from the Relativistic 
Heavy Ion Collider (RHIC) within the context of a sensitive, power-law representation 
introduced previously by Trainor and Prindle (TP). Uncertainties in the Glauber 
model input and in the mid-rapidity multiplicity frequency distribution data are 
reviewed and estimated using the TP centrality methodology. The resulting errors in 
model-dependent geometrical quantities used to characterize heavy ion collisions {i.e. 
impact parameter, number of nucleon participants Npart, number of binary interactions 
A^bin, and average number of binary collisions per incident participant nucleon are 
presented for minimum-bias Au-Au collisions at y/sj^ = 20, 62, 130 and 200 GeV 
and Cu-Cu collisions at y/sJfN = 62 and 200 GeV. Considerable improvement in the 
accuracy of collision geometry quantities is obtained compared to previous Monte Carlo 
Glauber model studies, confirming the TP conclusions. The present analysis provides 
a comprehensive list of the sources of uncertainty and the resulting errors in the above 
geometrical collision quantities as functions of centrality. The capability of energy 
deposition data from trigger detectors to enable further improvements in the accuracy 
of collision geometry quantities is also discussed. 
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1. Introduction 

Observables in relativistic heavy ion experiments are often reported as functions of 
measured total inelastic cross section fraction (centrality) and are related to the initial 
collision geometry using Monte Carlo Glauber (MCG) or optical Glauber models [HEIIS]. 
An example is charged particle multiplicity Nch versus number of participating nucleons 
Npart- It is therefore imperative that the level of accuracy of the Glauber models be 
understood and that relevant experimental information be used to estimate collision 
centrahty. For instance, previous studies [H [21 H] concluded that the uncertainties in 
relating peripheral collision data to the corresponding initial collision geometry were very 
large, thus contributing to the omission of those data from publications [1]. In another 
example, studies of high-p^ suppression of charged particle production for peripheral 
to mid-central collisions are presently limited by large uncertainties in the estimated 
number of binary collisions which is used as a scaling variable [5j. 

In this work we re-examine the accuracy of Monte Carlo Glauber model descriptions 
of minimum-bias multiplicity frequency distribution data from RHIC using updated 
density and multipliciity production model parameters and a sensitive power-law 
inspired representation introduced by Trainor and Prindle (TP) [6]. After demonstrating 
the accuracy of our Monte Carlo Glauber model in describing RHIC data we then 
use the power-law analysis method of Ref. [6] to estimate the uncertainties in the 
mapping relationships between initial stage collision geometry quantities and centrality 
for heavy ion collision systems relevant to the RHIC program. We find that previous 
uncertainties [H [21 SI in centrality measures are too pessimistic. 

The centrality analysis method of Trainor and Prindle [6] exploits the approximate 
power-law dependence of the multiplicity frequency distribution data dNevt/dNch {Nevt 
is the number of triggered events) for minimum- bias collisions and uses proton-proton 
(p-p) multiplicity production data to constrain the peripheral collision end-point of the 
dNevt/dNch distribution. Their analysis demonstrates that uncertainties in centrality 
determination and MCG geometry measures can be significantly reduced, particularly 
in the peripheral region. 

In this paper MCG results and errors for centrality bin average quantities (b) (mean 
impact parameter), (Npart), (A^bm), and u = {Nun) / {{Npart) /'^) [7] are presented for 
minimum-bias Au-Au collisions at ^/snn = 20, 62, 130 and 200 GeV and for Cu-Cu 
at 62 and 200 GeV. In this study centrality is based on charged particle production 
at midrapidity. Other centrality definitions appear in the literature but will not be 
considered here. However, the MCG model and analysis methods presented here and in 

X A minimum-bias trigger typically refers to the detection of forward going spectator fragments 
from both colliding nuclei plus a minimal requirement for particle production transverse to the 
beam direction. Each of the four RHIC experiments utilize a common lowest level trigger detector 
system based on calorimetric detection of neutrons at zero degree scattering angle using two zero- 
degree calorimeters (ZDC) placed symmetrically upstream and downstream from the beam-beam 
intersection region. The minimum-bias trigger systems for the STAR, PHENIX, PHOBOS and 
BRAHMS experiments are described in Nucl. Instrum. Meth. A 499, Nos. 2-3 (2003). 



Applicability of Monte Carlo Glauber models to relativistic heavy ion collision data 3 

Ref. [B| can be readily extended to any practical centrality determination method. 

The MCG model and updated input parameters are discussed in Sec. |2l The 
accuracy of the model for describing dNevt/dNch data from RHIC is demonstrated in 
Sec. 131 MCG predictions and errors are presented and discussed in Sec. HI Comparison 
of the MCG results and analytic parametrizations from Ref. [6] are presented in Sec. [5l 
In Sec. [6] we demonstrate via simulations how energy deposition data from transverse 
particle production, which are generally available at the trigger level from the RHIC 
experiments [8], can be used to further reduce the impact of background contamination 
and to mitigate the effects of event loss due to trigger and collision vertex finding 
inefficiencies. A summary and conclusions are presented in Sec. [71 Computational 
details are contained in three appendixes at the end. 

2. Monte Carlo Glauber Model 

The Monte Carlo Glauber collision model used here is based on a standard set of 
assumptions [3] which are appropriate for ultra-relativistic heavy ion collisions. These 
assumptions include the characterization of the collision in terms of a classical impact 
parameter (6), straight-line propagation of each incident nucleon through the oncoming 
nucleus, and a fixed transverse interaction range determined by the nucleon-nucleon (N- 
N) total inelastic cross section (ainei) in free space. The impact parameter was selected 
at random and the positions of the nucleons relative to the geometrical centers of the 
colliding nuclei were randomly distributed according to a spherical density p(r). Center- 
of-mass constraints were not imposed on the nucleon positions. Nucleon pairs in the 
colliding nuclei were assumed to interact hadronically if their relative impact parameter 
was < \/ainei/7r. 

Charged hadron multiplicity was assigned using the phenomenological two- 
component model of Kharzeev and Nardi p] for "soft" plus "hard" particle production 
processes where the mean number of charged hadrons in the acceptance (Nch) per unit 
pseudorapidity (77) was computed according to 



Parameters Upp = dNch{pp)/dr] and x depend on collision energy. Event-wise 
multiplicities were obtained by sampling the Gaussian distribution [9| 



where a is a multiplicity fiuctuation width parameter. Alternate distributions [e.g. 
Poisson, negative binomial produce quantitative effects on the multiplicity 

frequency distribution for very peripheral collisions but do not affect the present 
estimates of systematic errors. Nbm and Nch were both required to be > 1 in order 
for the simulated collision to be used in the analysis. Centrality bins were defined using 
the multiplicity frequency distribution dN^vtl dNch- Centrahty bin average quantities. 




(1) 
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(6), {Npart), {Nbin) and u = 2{Nhin) / {Npart) , Were calculated using the events within 
each bin. The acceptance for this study was |?7| < 0.5 and full 27r in azimuth. 

2.1. Matter Densities 

Monte Carlo Glauber simulations require the distribution of the centers of the nucleons 
in the nuclear ground state, Ppt,m{r), the point-matter density. For ^^Cu and ^^'^Au 
these were estimated using the measured charge densities [H] and the Hartree-Fock 
calculations of Negele using the Density Matrix Expansion (DME) framework [12] for 
the neutron - proton density differences. The charge densities for ^^Cu and ^^^Au were 
represented by a Woods-Saxon distribution, 

p{r) = po{l + exp[ir-c)/z]}-\ (3) 

where the radius and diffuseness parameters are listed in Table [H The point matter 

Table 1. Charge and point matter density Woods-Saxon parameters for ^'^Cu and 
-'^^^Au in fm. 



Density 
Parameter 


63 Cu 

EmpiricaP 


DME^ 


Empirical'' 


DMEi^ 


^chrg 


4.214±0.026 


4.232 


6.38±0.06 


6.443 


^chrg 


0.586±0.018 




0.535±0.027 




\ crirgl 


3.925±0.022 


3.899 


5.33±0.05 


5.423 


\' pt,mi 


4.195±0.085 
0.581±0.031 
3.901 


4.213 
3.875 


6.43±0.10 

0.568±0.047 

5.41 


6.495 
5.502 



^ Charge density results based on electron scattering analysis [TT] ; estimates of point 
matter densities as discussed in the text. 
^ Density Matrix Expansion predictions [T2]. 

densities assumed for the present analysis are spherically symmetric with a Woods-Saxon 
radial distribution where the half-density and rms radii were estimated by adding the 
Hartree-Fock DME point matter - charge distribution differences to the measured charge 
density radii where, 

/^2 \l/2 ^ / 2 \l/2, r/2 \l/2_/2 U/21 

v pt,ml \' chrgi ' \_\' pt,ml v chrgi J DME ' ^ ' 

where Cchrg and ('"c/irg)^^^ denote the measured radii from electron scattering while the 
quantities in the square brackets are the DME predictions. The diffusivity parameter 
Zpt,m was obtained from the quantities in Eq. (jlj) and the approximate relation, 
(r^) = |c^[l + |(vr2;/c)^] [l3]. The nominal radius and diffusivity parameters assumed 
here for Au and Cu are listed in Table [TJ 

The uncertainties in Cpt^m and Zpt^m were obtained by summing the independent 
errors in the proton [11] and neutron point matter densities in quadrature. The 
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latter is constrained by theoretical and experimental information about the neutron 
- proton density difference. Theoretical nuclear structure model predictions for the 
neutron— proton rms radii differences for isotopes in the Cu and Au mass range agree to 
within about ±0.06 fm [Hj. In general the theoretical predictions agree with analyses of 
medium energy proton-nucleus elastic scattering data [15] which are typically uncertain 
by about ±0.07 fm for isotopes in the Cu and Au mass range. The uncertainty in the 
neutron density rms radii relative to the proton density was therefore assumed to be 
±-\/0.06^ ± 0.07^ fm ^ ±0.09 fm and the corresponding uncertainty in the matter rms 
radii was ±{N/A)0.09 fm, where N,A are the number of neutrons and nucleons in the 
isotope. Theoretical contributions to the errors in Cpt^m and Zpt^m were conservatively 
estimated by requiring each to independently account for the ±(N/A)0.09 fm error. The 
latter theoretical errors were combined in quadrature with the corresponding errors for 
Cchrg and Zchrg from analysis of electron scattering data to obtain the final errors listed 
in Table [H 



2.2. N-N Inelastic Cross Section 

The N-N total inelastic cross sections used here were based on total cross section 
measurements for p-p collisions (±1 mb uncertainty) and elastic total cross sections 
for p-p and p-p (±0.5 mb error) [16]. Proton-neutron total cross section data are not 
available in the energy range studied here. The results for energies ^/s = 20, 62, 130 and 
200 GeV are respectively 33, 35.3, 38.7 and 41.7 mb, each being uncertain by ± 1.1 mb. 

It is possible that the effective interaction cross section between colliding nucleons 
inside a nucleus differs from that in free space (density dependence) or that the strength 
and range of the effective N-N interaction changes with each successive collision as in 
the "used" nucleon scenario [IT] or the "strict" participant scaling model [6] (limiting 
case of the used- nucleon model in which nucleons interact only once). The study of 
density dependent effects is well beyond the scope of the present analysis. 



2.3. Two-Component Multiplicity Production Model 

Parameters Upp = dN^hipp) / drj\ri=Q at y/s = 62, 130 and 200 GeV are 2.01, 2.25 and 
2.43 (±0.08 error for each), respectively, using the energy dependent parametrization 
of the UA5 [H] and CDF data given in Ref. [19]. For the ^ = 20 GeV data, which is 
outside the energy range parametrized in [19], the average of ISR [20] and FNAL [21] 
measurements summarized in Ref. [22] was assumed where Upp = 1.4 ± 0.12. 

Binary scattering parameter x in Eq. ([T]) was estimated by fitting {dNch/ drj) / {O.bNpart) 
data from STAR [1 123], PHENIX [21], and PHOBOS [25l |26] versus z/ with Eq. (P 
rewritten as, 

A^^ = "-[' + "^'^-')]' 
assuming the above values for ripp. The resulting values of x are 0.07, 0.09, 0.09 
and 0.13 (±0.03 errors for each) for the 20, 62, 130 and 200 GeV data, respectively. 
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Analysis of the 19.6 GeV Au-Au data by the PHOBOS experiment assumed 
Upp = 1.27 ± 0.13 [20] and obtained x = 0.12, leading to claims in the literature that 
hard scattering contributions to multiplicity do not change with collision energies from 
20 to 200 GeV [251 ISSl [271 [2H] • A more recent compilation [22| of Upp measurements in 
this energy range indicates a larger value of Upp = 1.4 ± 0.12. Fitting Eq. ([5]) to the 
combined PHENIX [21] and PHOBOS [25] 19.6 GeV data with Upp = 1.4 resulted in a 
smaller x = 0.07. PHOBOS 62 GeV Au-Au data [26], when plotted versus z/, do not 
linearly extrapolate to Upp in contrast to what is expected from Fig. 4 of Ref. [26] and 
Fig. 1 of Ref. [28] from the PHOBOS experiment. However, 62 GeV Au-Au results from 
the STAR experiment linearly extrapolate to Upp at z/ = 1, resulting in the fitted value 
X = 0.09 used here, x parameters from the three RHIC experiments generally agree for 
the 130 and 200 GeV data; the PHOBOS values of x = 0.09 [25] and 0.13 [25] at 130 
and 200 GeV, respectively, were confirmed and used here. It is possible that parameter 
X is affected by other processes in addition to hard and semi-hard partonic scattering 
in the initial collision stage where the latter mechanisms would be expected to follow a 
log y/s dependence. 

The variance of the multiplicity distribution V{Nch, Nch) in Eq. ([2]) is given by 
aNch where values of parameter a < 1 or > 1 represent multiplicity fluctuation 
suppression or excess, respectively, relative to pure statistical fluctuations (a = 1). 
In general a non-vanishing integral [29l [30l [3T] of two-particle correlations [32] over 
the acceptance requires a ^ 1. In principle, correlation measures of the type 
reported in Ref. [32] could be used to determine parameter a using the relationship 
between fluctuations and correlations developed in Refs. [29l [31] and discussed in 
Appendix A. On the other hand, the Kharzeev and Nardi two-component multiplicity 
model with distribution V{Nch, Nch) constitutes a phenomenology for describing event- 
wise multiplicity frequency distribution data. The phenomenological approach based on 
fits to multiplicity distribution data, shown in the next section, was used to estimate 
parameter a in the present analysis. The uncertainty in the width of V{Nch, N^h) was 
estimated by fitting the data. The uncertainty in the analytic form of the multiplicity 
fluctuation distribution was accounted for by comparing Monte Carlo Glauber results 
assuming the Gaussian V {Nch, Nch) in Eq. ([2]) (with a = 1) with results assuming a 
negative binomial distribution (NBD) [10] in place of Eq. ([2]) as explained in Appendix A. 

3. Fits to 130 GeV Au-Au Data 

The 130 GeV Au-Au minimum-bias negatively charged hadron multiplicity frequency 
distribution data from STAR [4] for 60K events were fit by adjusting parameters Upp (for 
negative hadrons) and a where x was fixed to 0.09 as discussed above. The acceptance 
was defined by transverse momentum (pt) > 0.1 GeV/c, \ri\ < 0.5, and A0 = 2n. 
The fit (solid histogram) is shown in the left panel of Fig. [1] in comparison with data 
(solid dots), where the optimum values of Upp and a are 1.110 ± 0.004 (consistent with 
Upp = 2.25 ± 0.08 for charged particle yields) and 1.04 ± 0.10, respectively. Similar fits 
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Figure 1. Negative liadron multiplicity frequency distributions for Au-Au minimum- 
bias collisions at y^J/vw — 130 GeV normalized to the total inelastic cross section in 
barns. Left panel (a): semi- log plot comparing STAR data [4] (solid dots) with the 
Monte Carlo Glauber model fit (histogram). Middle panel (b): same data and Monte 
Carlo fit on a log-log plot showing the power-law dependence and the exponent (slope) 
of approximately -3/4. Right panel (c): same data (solid dots) and Monte Carlo fit 
(solid line) plotted as da/dN^['^ versus N^-'^. 



(not shown) to the charged particle minimum-bias muhiphcity distribution from the 
same data set resulted in a = 0.94 ± 0.15. Parameter a was therefore set to 1.0 ± 0.2 for 
all four energies. 

In the TP analysis [6] it was shown that minimum-bias multiplicity frequency 
distributions for relativistic heavy ion collision experiments and Monte Carlo Glauber 
models approximately follow a power-law distribution. This is illustrated by plotting the 
data and MCG fit from the left panel of Fig.[T]on log- log axes in the middle panel. Except 
near the end-points the data can be described to within 10% with a slope (exponent) of 
approximately —3/4. The collision event yield is approximately proportional to N^^^^. 
Therefore distribution da/dN^^ versus N^_^ is approximately constant, where 

da dNh- da 3/4 da . 
777 = VTJ-, = ^N' K, const. (6) 

dNlL' dNl'.' dN,- ^ dN,- ^ > 

as shown in the right-hand panel of Fig. [1] for data (solid dots) and MCG fit (solid 
line). The power-law plotting format in the right-hand panel enables a more sensitive 
comparison between the model fit and data than in the usual semi-log format (left 
panel). The MCG model fit is consistent with the data except for the first two data 
points at low N^^- ■ Except near the end-points the da/dN^L^ data are constant on N^[^ 
to within an overall variation of 20% in comparison to the da/dN^- data which span 
nearly two orders of magnitude within the same multiplicity range. 

The results in Fig. [1] demonstrate the efficacy of the Monte Carlo Glauber model 
with two-component multiplicity production for accurately describing the measured 
multiplicity frequency distributions at RHIC. Based on this outcome we conclude that 
the present model is reasonable to use in estimating centrality bin average quantities 
and their systematic errors for RHIC data. 
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4. Results and Error Analysis 

Ensembles of one- million, minimum-bias {i.e. random impact parameter) Monte Carlo 
collisions were generated for each of the six systems studied here. Centrality bin averaged 
quantities (Nch), {b), {Npart), {Nbm) and z/ are hsted in Tables [2H7] using the nominal 
parameter values discussed above. Statistical errors are typically 0.1-0.2% and always 
< 0.5% of the nominal bin averages and in all instances are much less than the systematic 
errors discussed below. Results in Table [2] for {Npart) for Au-Au collisions at 200 GeV 
are systematically larger than published Monte Carlo Glauber predictions in Ref. [1] by 
approximately 5% for peripheral centrality bins. This systematic increase is primarily 
caused by multiplicity fluctuations in the present model and basing centrality on Nch 
here rather than on Npart as was done in Ref. [1]. 

Average charged particle multiplicity per participant pair for |?7| < 0.5 at mid- 
pseudorapidity as a function of centrality {u) is shown in Fig. [2] for Monte Carlo 
simulated Au-Au colhsions at ^/snn = 200 GeV (solid dots). The p-p limit, Upp = 
2.43±0.08, is indicated by the solid square symbol. The data display a linear dependence 
on (z^ — 1) as expected from Eq. ([5]) except for the most-peripheral 90-100% centrality bin 
where multiplicity fluctuations significantly reduce the average Nch- The slope agrees 
with UppX = 0.32; the linear extrapolation to the p-p limit is evident. 



CM 



n 
a. 



Au-Au 200 GeV 
p-p 200 GeV 



Figure 2. (color online) Monte Carlo Glauber results for {Nch) / {Npart l"^) for one- 
million Au-Au collisions at ^/smn — 200 GeV using the nominal parameters discussed 
in the text (solid dots). Linear fit (dashed line) to the Au-Au results (excluding the 90- 
100% centrality bin) accurately extrapolates to the p-p multiplicity from UA5 [TSlllQj . 
Monte Carlo data for Au-Au include statistical errors only which are smaller than the 
symbols. 
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4.1. Error Estimation Method 

The sources of uncertainty which cause systematic error in the colhsion geometry 
quantities can be organized into three categories. The first includes fitting uncertainties 
that occur when the corrected, multiphcity frequency distribution data are described 
by adjusting the parameters of the multiphcity production model in Eqs. ([I]l2l). The 
second includes the uncertainties in the matter density, nucleon-nucleon total inelastic 
cross section, and functional representation of the multiplicity fluctuation distribution 
V{Nch, Nch)- The third includes the uncertainty in the corrected multiplicity frequency 
distribution data dNf^^t/ dNch- The latter arise from trigger inefficiency, collision 
vertex finding inefficiency, background contamination, and particle tajectory finding 
inefficiency. 

Systematic errors in centrality bin average quantities due to fitting uncertainty were 
estimated by comparing the nominal results to that obtained with parameters x and 
a in Eqs. ([lH2]) individually varied within their respective fitting errors (Sees. 12.31 and 
|3]). Uncertainty in parameter npp directly affects {Nch) but does not affect (6), {Npart), 
{Nhin) or u. 

Error estimates due to uncertainties in the matter density radius and diffuseness 
parameters and the N-N inelastic cross section require a refitting of the dNevt/dNch 
distribution. This was accomplished via a phenomenological adjustment of the 
multiplicity production model described in Appendix B. Fit recovery via |xP 
minimization for MCG simulations with ample statistics (of order 10^ collisions) is 
computationally intensive. The method in Appendix B is fast and accurate. Errors due 
to the uncertainty in the mathematical representation of V{Nch, N^h) were estimated 
by comparing the nominal MCG results using the Gaussian distribution in Eq. ([2]) with 
results assuming a negative binomial distribution [181 [IH] as explained in Appendix A. 

Uncertainties in the dNevt/dNch data arising from trigger inefficiency, collision 
vertex finding inefficiency, and background contamination mainly affect the low 
multiplicity region. Minimum-bias trigger efficiencies at RHIC are 92.2l3Q% [211 [33] , 
94 ± 2% [4J, 96% [H] and 97 ± 3% [25l [351 [36] . Trigger inefficiency causes the lower 
Nch half- max point of the dNc^t/dN^l'^ distribution to shift to larger N^l'^. The position 
of this point for the dNevt/dNH^ distribution when corrected for trigger inefficiency 
(assuming the power-law behavior) has a relative uncertainty of about ±10%. [§| However, 
knowledge that the lower half-max point is constrained by p-p scattering allows the 
relative uncertainties of the lower half-max point on N^^ to be reduced to 1/4 of the 
uncertainties in Upp, or to ±0.8%, ±0.9%, ±1%, and ±2% for the 200, 130, 62 and 
20 GeV data, respectively. 

Primary collision vertex reconstruction efficiency is approximately 100% for events 
with Nch of order a few tens and greater but falls precipitously for Nch < 10 [37] . 

§ Trigger efSciencies are uncertain by about ±2 to ±3%. The fractional uncertainty in the lower half- 
max position is ±(0.02 to 0.03)(iVi/^„,, - A^'/.tuJ/^cY,™™ « 10% where iV^/t.™ and iV.Y,!,, are 
the lower and upper half-max points of the dN^yt/ dN^j^^ distribution. 
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PHOBOS [251 and STAR [U [37] both report uncertainties in their overall vertex 
finding efficiencies of about ±2% for minimum-bias collisions. Vertex finding inefficiency 
increases the slope and half- max position of the low multiplicity edge of the dNevt/dN^l'^ 
distribution. An uncertainty of ±2% in the overall efficiency results in about ±10% 
uncertainty in the low N^l'^ half-max position. However, the p-p data constrain this 
uncertainty. In this analysis we assumed 100% vertex finding efficiency for collisions 
producing N^h > 14 [37], or N^l^ > 1.93, and allowed the slope of the lower edge of 
the dNevt/dN^[^ distribution to vary such that the half-max position varied by ±0.8%, 
±0.9%, ±1%, and ±2% for the 200, 130, 62 and 20 GeV data, respectively, as in the 
preceding paragraph. 

The principle sources of background contamination are from ultra-peripheral two- 
photon interactions |38] and beam-gas collisions. The former process corresponds to 
coherent photon-photon interactions which excite both nuclei, followed by neutron decay 
(which activates the minimum-bias trigger detectors) and accompanied by resonance (s) 
production which decays into charged particles transverse to the beam direction. 
Transverse particle multiplicities from UPC events are typically < 2 {e.g. p-meson 
decay) and generally < 4 [38] for \rj\ < 1 at midrapidity. UPC backgrounds are 
therefore restricted to the most-peripheral (90-100%) centrality bin. UPC yields should 
be approximately proportional to (^1^2)^ [38J (charge numbers for colliding ions 1 and 
2) whereas beam-gas contamination should scale with beam current. Other background 
events, e.g. mutual Coulomb dissociation processes such as7±A^ A* ^B±n 
for both nuclei, can be eliminated by requiring minimum transverse particle production. 
Remaining background events will appear near the lower Nch edge of the dNg^t/dN^l^ 
distribution. 

Estimates of background contamination in minimum-bias trigger data at RHIC 
range from 1 ± 1% [33] to 6% lU H] overall for Au-Au at ^/sNN =130 GeV corresponding 
to 0-20% and 60%, respectively, of the hadronic collision event yield in the 90-100% 
centrality bin. Most of the UPC events occur at N^h < Upp and can be eliminated 
by cuts on the number of transverse charged particles. The remaining background 
contributions for N^h > npp are less than the amounts listed above. For the present 
analysis background contamination was assumed to be dominated by UPC events and 
to diminish in magnitude with collision energy and (^1^2)^. In the present analysis 
background contamination was applied to the nominal 90-100% centrality bin assuming 
3%, 2% and 1% overall contamination levels in the Au-Au data at ^/SNN = 200, 130 
and 62 GeV, respectively. UPC contamination for Au-Au collisions at 20 GeV and 
for Cu-Cu at 200 and 62 GeV was estimated to be negligible. However, calculations 
were done for the latter three systems assuming a 1% overall background contamination 
(10% contamination within the nominal 90-100% centrality bin) in order to provide a 
reference for further systematic error estimation. 

The above sources of systematic uncertainty primarily affect the position, slope 
and shape of the dN^vt/ dN^^^ distribution at the lower end-point. These changes 
impact the MCG model which must describe those data and, in turn, the collision 
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geometry quantities like {Npart)- The phenomenological method in Appendix B was 
used to estimate the changes in the colhsion geometry quantities relative to the nominal 
values. The latter differences were taken as the estimated errors. 

Background levels in collider experiments vary significantly depending on beam 
quality, beam current and interaction rate. Excessive trigger backgrounds beyond that 
considered here may result from beam-gas interactions during periods of high integrated 
beam currents. Collision event pile-up in the particle tracking detectors during periods 
of high luminosity adversely affect collision vertex finding. Either condition may be 
so severe as to preclude access to the low multiplicity range of the minimum-bias 
distribution. Even so, the power-law dependence and p-p end-point constraints enable 
accurate centrality estimates and collision geometry assignments to be made for the 
remaining minimum-bias data. 

Charged particle trajectory reconstruction efficiencies in the large acceptance 
tracking detectors at RHIC are typically 70 - 95% [U [2ll [25] and decrease approximately 
linearly with particle density in the detectors by about 20% from most-peripheral 
to most-central collisions p57]. Uncertainties in the assumed track reconstruction 
efficiencies were estimated by comparing corrected data for {dNch/ drj\r^=o) / {Npartl'^) 
versus Npart between the RHIC experiments for Au-Au collisions at 20 GeV [211 [25] 
130 GeV [H, [211 [25] and 200 GeV [211 [25]. The comparisons indicate uncertainties in 
both the overall tracking efficiencies and in the dependence on particle density in the 
tracking detectors. The latter variation in efficiency from peripheral to central collisions 
is about 20 ± 8%. 

Overall changes in tracking efficiency are compensated in the MCG model by 
multiplicative adjustments to parameter Upp and therefore have no effect on the 
centrality measures reported here. Changes in the assumed tracking efficiency 
dependence on N^^ affect the lower and upper end-point positions of the dN^vt/ dN^l^ 
distribution and distort its shape. The distortions must be accounted for by the MCG 
model in order to determine the net effect on the collision geometry measures. The effects 
of the 8% uncertainty in the multiplicity dependence of the trajectory reconstruction 
efficiency were estimated by generalizing parameter Upp in Eq. ([1]) to ripp(l + aNpart) 
where a = ±0.00023 and aNpart = ±0.08 for most-central Au-Au collisions. The same 
value for a was assumed for Cu-Cu. Systematic errors were estimated as the differences 
between the nominal centrality bin averages and those obtained assuming a = ±0.00023. 
The choice to maximize multiplicity variation for most-central collisions is arbitrary. As 
a result systematic errors in centrality bin average multiplicities were not included in 
Tables [2}[7l However, this ambiguity does not affect the resulting systematic errors in 
the other centrality measures reported here. 

4-2. Error Results 

The combined systematic errors (all components added in quadrature) in both 
magnitude and relative percent (given in parentheses) are listed in Tables [2] - [3 for the 
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six collision systems studied here. Impact parameter uncertainty is about ±(2 to 3)% 
for Au-Au and ±(2 to 6)% for Cu-Cu. Uncertainty in (Npart) is about ±(4 to 8)% for 
Au-Au and ±(3 to 7)% for Cu-Cu for peripheral centralities, reducing to about ±1% or 
less for central collisions where full geometrical overlap of the two nuclei suppresses the 
dependence of Npart on variations in the nuclear surface geometry. Systematic errors 
in {Nhin) vary from about ±(4 to 12)% for Au-Au colhsions and ±(5 to 10)% for Cu-Cu 
for central to peripheral collisions. Errors in u (±4% or less) are suppressed due to 
covariation of (Npart) and {Nbm). 

Individual contributions to the systematic errors in collision geometry bin averages 
(6), {Npart), {Nbin) and z/ for Au-Au and Cu-Cu colhsions at ^/snn = 200 GeV are 
summarized in Tables [8] and [9l respectively. Similar results were obtained for the other 
four collision systems. Errors are given in percent where the absolute values of positive 
and negative errors were averaged together. The three numbers listed for each instance 
correspond to the average errors within the 60-100%, 20-60% and 0-20% centrality bins, 
respectively. The dominant errors are due to uncertainties in the matter density and 
the N-N total inelastic cross section. Errors due to uncertainties in the analytic form of 
the phenomenological multiplicity fluctuation model are only significant for peripheral 
collisions. Errors due to uncertainties in the multiplicity dependent particle track 
reconstruction efficiency and in the trigger and vertex reconstruction inefficiencies are 
negligible when constrained by p-p data as shown previously by Trainor and Prindle [6]. 

Percent uncertainties due to possible background contamination are listed for all six 
collision systems in Table [THl The largest errors occur in the most-peripheral centrality 
bin as expected. The absolute magnitudes of the errors decline rapidly with increasing 
centrality and are negligible for the centrality bins not listed in the table. Reference 
errors for Au-Au at 20 GeV and Cu-Cu at a/J/vjv = 200 and 62 GeV were based on an 
assumed 10% background contamination in the 90-100% centrality bin. 

Estimates of total systematic error when background contamination differs from 
that assumed here can be obtained by removing the error contributions in Table [10] 
from the total errors in Tables [2] - H] for Au-Au collisions at 200, 130 and 62 GeV, 
respectively, and then adding (in quadrature) the appropriately scaled background 
errors from Table [TOl For Au-Au collisions at 20 GeV and Cu-Cu collisions at 200 
and 62 GeV the scaled errors from Table [10] should be combined in quadrature with the 
total systematic errors in Tables [5] - [7] If much larger backgrounds than those assumed 
here are encountered, then the present Monte Carlo results should not be scaled; rather 
the errors should be recalculated. 

Overall, the systematic errors in collision geometry bin averages for non-peripheral 
centralities are dominated by uncertainties in the nuclear geometry and o"j„e«- Errors in 
the more peripheral bins are dominated by background contamination and ambiguities 
in the analytic form of the multiplicity fluctuation model. 
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5. Analytic Parametrizations 



The power-law description of heavy-ion collision centrality developed by Trainor and 
Prindle |6] prescribes simple, analytic parametrizations of the N parti Nhm, ^ and Nch 
dependences on total inelastic cross section fraction a joQ. For N^art the running integral 
relation on (1 — a/ao) is accurately given by [6] 

[^~'^/'^o) — i 11 [i ) 

{^■^paTt,max / '^') i_-^part,min/ '^^ ^ 



where 



Similarly 



(A/part/2)4 = {Npart,min/2)^a/ao + (A/part,max/2) * (1 - (j/ao). (8) 



where subscripts min and max refer respectively to the lower and upper half-max end- 
point positions of the da/d{Npart/'^y^^ and da/dN^-^ distributions. 

The above parametrizations are shown as the solid lines in the two upper panels 
of Fig. [3] in comparison with the present MCG model results (solid dots) from Table [2] 
for Au-Au collisions at ^/sNN = 200 GeV. The end-point parameters for the power- 
law parametrizations were Npart,min/'2 = 0.75, Npart,max/'2 = 189, Nun,min = 0.75 and 
Nbin,max = 1144. The parametrization for u = Nbm/ {Npartf^) is compared with the 
MCG result in the lower- left panel. The average multiplicities using Eq. ([5]), the values 
for Upp and x from Sec. 12. 3[ and the above parametrizations for Npart and u (solid 
curve) is compared with the MCG result in the lower-right panel. The simple power- 
law parametrizations introduced in Ref. [6] quantitatively describe the MCG results, 
thus confirming their utility and precision. 

Calculation of the analytic power-law parametrizations require the upper half-max 
end-point values for Npart and N^in from the MCG model. These quantities are listed 
in Table [TT] for the six collision systems studied here. 



6. Discussion 



Very peripheral collision data from the four RHIC experiments have generally remained 
unpublished due to concerns with possibly large, and not well understood background 
contamination, trigger inefficiencies, and primary collision vertex finding inefficiencies. 
Trainor and Prindle [6] originally showed, and the present analysis confirms, that the 
power-law dependence of the multiplicity frequency distribution data together with 
knowledge of the proton-proton multiplicity enables accurate centrality information to 
be obtained in spite of these uncertainties. Nevertheless, UPC events in principle cause 
the lower end-point portion of the multiplicity frequency distribution for minimum-bias 
A-A collisions to differ from the p-p limit and these backgrounds will contaminate the 
A-A data in the most-peripheral centrality bins. In this section we discuss additional 
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Figure 3. Comparison of Monte Carlo Glauber centrality quantities (solid dots) 
for Au-Au minimum-bias collisions at ^ysNN = 200 GcV from Table [5] with analytic 
parametrizations from Trainor and Prindle [5] based on the power-law description of 
heavy- ion collisions as discussed in the text. The four panels show from upper-left to 
lower-right the dependences of N^/n^ v and on relative cross section 

fraction (1 — a jao). 

analysis methods using data available to the RHIC experiments to minimize backgrounds 
and to provide contamination level estimates for the accepted centrality bins. 

Transverse particle production information via scintillators, silicon detectors, 
and/or calorimetry are available in the RHIC experiments at the trigger level [33 SHI 
132] . These data record total energy deposition in the sensitive detector material for 
each (minimum bias) triggered collision event including contributions from A-A hadronic 
collisions, ultra-peripheral coUisions, beam-gas and other backgrounds, etc. but do not 
include particle track finding and primary collision vertex finding inefficiencies. If the 
integrated yields from these detectors, denoted as SUM, is approximately proportional 
to Nch, then the frequency distribution will follow an approximate power-law such that 
the data can be usefully represented by dNtrig/dSUM^^^ versus SUM^^^, where Ntrig is 
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the number of triggers. 

MCG simulations for the integrated yields from the STAR central trigger barrel 
(CTB) [39] plastic scintillator detector {\rj\ < 1, 27i azimuth coverage) were done for p-p 
and Au-Au minimum-bias collisions at ^/SNN = 200 GeV. Details of this simulation are 
discussed in Appendix C. The results for dNtrig/ dSUM^^^ versus SUM^^^ are shown in 
Fig. m Both the p-p and Au-Au simulations reproduce the general shapes of measured 
STAR CTB minimum-bias trigger yields [8]. The simulated Au-Au CTB yields follow 
the SUM~^/^ power-law distribution very well where the lower half-max point coincides 
with the mode of the p-p distribution. Trigger inefficiency in Au-Au causes a depletion 
at low multiplicity as indicated by the green, dashed curve. Over- or under-corrected 
yields result in the blue dotted curves. From Fig. H] it is clear how the p-p trigger yield 
can be used to constrain corrections for trigger inefficiency in A-A. 




Figure 4. (color online) Monte Carlo simulation results for the minimum-bias trigger 
frequency distribution for transverse particle production integrated trigger detector 
yield quantity SUM to the 1/4 power. Simulated yields are shown for IM 200 GeV 
Au-Au (solid black curve, multiplied by 3) and IM proton-proton (solid red curve) 
collisions. The yield for Au-Au collisions when trigger inefficiency is large is illustrated 
by the lowest (right-most) green dashed curve (hand-drawn sketch). Over- and under- 
corrected yields are similarly illustrated by the upper (left-most) and middle dotted 
blue lines, respectively. Background contamination contributes at lower multiplicity 
as illustrated by the dashed-dotted magenta curve (hand-drawn sketch). 

Background triggers appear at lower values of SUM where UPC events typically 
produce of order two charged particles at mid-rapidity compared with the average from 
p-p collisions of about 5 in jr^l < 1 and therefore appear as a peak or enhancement 
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(dashed-dotted curve in Fig. H]) below the mode of the p-p distribution. Comparisons 
between the measured p-p and Au-Au dNtng/dSUM^^^ distributions and the power-law 
simulations for Au-Au provide a reasonable means for defining additional event cuts 
to reduce background contamination. This information can also be used to estimate 
the contamination level for accepted events. Measurements of the dependence of the 
lower end-point region of the dNt^ig/ dSUM^^^ distribution on integrated beam current, 
luminosity, and ion species would help disentangle background contributions from beam- 
gas collisions, event pile-up, and UPC events. A similar power-law analysis using 
reconstructed particle tracks, but without a primary collision vertex requirement, would 
permit the vertex finding inefficiencies to be estimated and corrected. 

7. Summary and Conclusions 

A Monte Carlo Glauber and two-component multiplicity production model with 
fluctuations for high energy heavy ion collisions was used to describe minimum-bias 
multiplicity frequency distribution data from RHIC. Updated Woods-Saxon radii and 
diffusivities for the point matter densities of ^^^Au and ^'^Cu were determined using a 
combination of charge density measurements from electron scattering and theoretical 
Hartree-Fock predictions. The binary scattering parameter x was estimated using a 
compilation of available RHIC data from 20 to 200 GeV. The values for x obtained here 
display some energy dependence, although with large uncertainty. The model was fltted 
to the 130 GeV Au-Au minimum-bias data [1] in both the conventional semi-log format 
and in the more sensitive power-law inspired format introduced in Ref. |6] . Quantitative 
agreement between the present model and data was obtained. 

Systematic errors in centrality bin averages for the geometrical quantities (5), 
(Npart), {Nbin) sjid V = {Nbin) / {{Npart) / 2) Were estimated based on charged particle 
multiplicity for |?7| < 0.5 and full 27r azimuth acceptance. The sources of systematic 
error considered here included uncertainties associated with the model (nuclear density, 
nucleon-nucleon inelastic total cross section, multiplicity production model analytic 
form and parameters) and with the corrected minimum-bias multiplicity frequency 
distribution data (background contamination, trigger and collision vertex flnding 
inefficiencies, and particle trajectory reconstruction inefficiencies). The TP centrality 
analysis method was applied to the MCG predictions and error analysis for six collision 
systems relevant to the RHIC heavy ion program. 

The analysis showed that significantly reduced errors in the collision geometry 
quantities result when the uncertainties in the minimum-bias multiplicity frequency 
distribution data are constrained by the empirical power-law behavior and the minimum- 
bias p-p collision data. The reduction in errors is in agreement with the original 
TP power-law analysis; the resulting errors are in general smaller than previous 
estimates [H El 0] which did not utilize the power-law and p-p constraints. The accuracy 
of simple parametrizations of centrality dependent quantities developed by Trainor and 
Prindle [6] was confirmed. We also discussed how particle production data at the trigger 
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level {e.g. scintillator hits and calorimeter energy depositions) can be used within the 
power-law and p-p constraint methodology [6] to define additional event cuts which 
minimize background contamination and which enable the detrimental effects of trigger 
inefficiencies to be reduced. 

The systematic errors presented in this paper demonstrate the accuracy of collision 
geometry quantities which can be achieved with RHIC data by exploiting the power-law 
behavior of the A-A data and using constraints from minimum-bias p-p collision data. 
The results presented here assume that backgrounds from UPC events, pileup, beam-gas 
collisions, etc. are not too large and that the minimum-bias trigger and collision vertex 
reconstruction efficiencies are not too small for peripheral collisions. It is intended that 
the MCG results presented here will serve as a useful resource for the RHIC community 
and that the analysis method developed by Trainor and Prindle [6] and applied in 
this paper, together with the trigger data analysis method discussed here, will enable 
accurate description of, and better access to the heavy ion peripheral collision data from 
RHIC. 
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Appendix A. 

Multiplicity frequency distributions for high energy minimum-bias proton-proton 
collisions are well described by a negative binomial distribution (NBD) [TF, T8j given by 

In Eq. OA. II) n and (n) are random and mean multiplicities, the variance is (^)(1 + 
{n)/k), the variance excess relative to Poisson statistics is {n)'^/k, and the multiplicity 
fiuctuation variance excess per final-state particle is {n)/k. The latter expression serves 
as an operational definition of parameter k. For p-p collisions both (n) and 1/k increase 
approximately linearly with ln(s) [18] as expected from pQCD cross sections for parton 
scattering. 

Multiplicity fiuctuation variance excess per final-state particle within a given 
acceptance is equal to the integral of two-particle correlations on relative momentum 
coordinates within the same acceptance [29l[30l[3T]. For Au-Au collisions at 130 GeV [32] 
and at 62 and 200 GeV [43] the integral of charged particle correlations for |?7| < 0.5 
and 2iT azimuth acceptance is dominated by a two-dimensional (2D) correlation peak 
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Table 2. Centrality-bin averaged collision geometry quantities and errors for Monte 
Carlo Glauber Au-Au collisions at y/s^^ — 200 GeV using the nominal model 
parameters discussed in the text. Centrality was based on Nch for |??| < 0.5 and 
full 2tt azimuthal acceptance. Estimated positive and negative systematic errors are 
listed as magnitudes and percentages (in parentheses) of the mean values. Errors in 
mean multiplicities were not computed (see text). 
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at small relative opening angles. The magnitude of this correlation integral increases 
roughly with [y — l)lnyj/v7v- Multiplicity fluctuation variance excess per flnal-state 
particle was therefore approximated by 



(n) (n) p ^ 1 



[«:-l)(z/-l) 



(A.2) 



(z/o_5% - 1) 

where quantity p — p] is the variance excess for p-p collisions, z/o_5% is the 

average number of binary collisions per incident participant nucleon for the most-central 
(0-5%) A-A data, and k is taken to be the ratio of the 2D correlation peak amphtude 
for most-central A-A collisions to that for p-p collisions. 

Values of {n)/k for p-p collisions were obtained from UA5 [ini [IE] using either 
measured values of k or the UA5 energy dependent flt function for k. (n) for |?7| < 0.5 
is given by the values of Upp assumed in this analysis. Simulations demonstrated that 
the values of k obtained by UA5 with very large rj acceptance remained approximately 
constant at the smaller rj acceptance assumed in this analysis. Values of quantities 
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Table 3. Same as Table [2] except for Au-Au collisions at ^s^vAf ~ 130 GeV. 
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(n)/A;[yi,p-p] were 0.09, 0.24, 0.40 and 0.51 for = 20, 62, 130 and 200 GeV, 
respectively, k for Au-Au collisions at 20, 62, 130 and 200 GeV were estimated to be 
4.8, 5.0, 5.4 and 5.5, respectively, from analysis of correlation data at 130 GeV [32] 
assuming Iua/s scaling. The remaining values for Cu-Cu at 62 and 200 GeV were 
respectively 3.4 and 3.7 assuming [y — 1) scaling. Finally, for each simulated collision 
with event-wise v and N^h from Eq. ([I]) NBD distribution Eq. (IATD was sampled to 
obtain the event-wise multiplicity. 

Appendix B. 

Fit recovery of the nominal dNevt/dNch distribution following shifts in the nuclear 
density parameters and aimi was achieved by constructing a perturbative correction, 
S{dNch/dri), which was added to the right-hand side of Eq. ([T]). Parameter shifts affect 
the collision distributions on Npart and N^in and in turn the dNevt/dNch and dNevt/dNch 
distributions. The fit recovery term is defined such that the nominal dNcvt/dNch 
distribution is maintained when the frequency distribution on Npart changes. Running 
integrals of the distributions dNevt/dNch and dNevt/ d{Npart/'^) produce a one-to-one 
correspondence between Nch and Npart /'^ as a function of the number of summed events. 
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Table 4. Same as Table [2] except for Au-Au collisions at ^s^vat = 62 GeV. 
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1.75; 


f0.05(3.0) 
-0.04(2.2) 


60-70 


29.0 


11 gi+0-21(1.7) 
XX. 24(2.0) 


26 6+'-'^^-'^ 

-^"•"-0.7(2.7) 


97 c-+2.2(8.1) 
^'•^-1.5(5.5) 


2.07; 


f0.07(3.6) 
-0.06(3.0) 


50-60 


52.3 


10 94+°-''^'^ 

J-'J.S''±_Q 22(2.0) 


9+1.5(3.3) 
^"-^-1.1(2.3) 


dU.O_3 2(5 6) 


2.46; 


f0.10(4.1) 
-0.09(3.5) 


40-50 


87.1 


q 00+0.17(1.7) 
^•""-0.18(1.8) 


7/1 9+1-9(2.6) 
'^•^-1.1(1.5) 


^'J^-8_5 3(4,9) 


2.90; 


f0.12(4.3) 
-0.10(3.6) 


30-40 


136.9 


o 71 +0.15(1.7) 
' -"--0.17(1. 9) 


1 1 9 /I +2-0(1.8) 
-"--"-^•^-1.3(1.2) 


190 2+^^-^(^-°) 

-"-^^•^-8.7(4.6) 


00+0.14(4.2) 
'^•"^°-0.12(3.7) 


20-30 


205.5 


7 o(-+0. 12(1.6) 
' "^"-0.15(2.0) 


162 6+^-°^^-^^ 

-"-"^•"-0.9(0.6) 


01 p- 9+16.6(5.3) 
OXj.z_-^2.3(3.9) 


3.88; 


f0.16(4.1) 
-0.14(3.5) 


10-20 


300.7 


. ^0+0.09(1.6) 
"-'•^"-'-0.12(2.1) 


990 9+1-9(0.8) 
-^^^•^-1.0(0.4) 


c-qo ^+23.6(4.7) 
dUO.U_2o.4(4,o) 


4.39; 


f0.17(3.9) 
-0.17(3.8) 


5-10 


395.5 


A qr, + 0.06(1.4) 

^•"^^-0.07(1.8) 


9QQ 0+1-3(0.5) 
zyo.o_Q 2) 


700 c-+30.6(4.3) 
(UO.d_2g 3(4 0) 


4.8O; 


f0.19(3.9) 
-0.18(3.8) 


0-5 


480.8 


9 r,o+0.05(2.4) 
^■'^°-0.06(2.4) 


0^d.U_2.o(0.6) 


88 4 9+38-8(4-4) 
°°^-^-39.5(4.5) 


5.13; 


f0.20(4.0) 
-0.21(4.0) 



Using the power-law representation this running integral relation is given by, 

dN'J''-^^ 4iV„,/2)'^ ^,f*°',„ (B.l) 



which defines a locus of ordered pairs [(A'part/2)^'''*, iV^/^], which can be compactly 
expressed discrete function defined by 

NIL" = ^ll%'^D[{Npart/2f'% (B.2) 
In Eqs. (IB.ip and (1B.2|) subscripts D and D' refer to the nominal (reference) distribution 
when D = R and to the shifted distribution when D = S. Lower limit Ncl'lnin — "^vv 
for the reference distribution. In order to recover the nominal dN^vt/ dNch distribution 
the correction term must be computed by the following difference, 

K^)=KKS^n[{^)^-KKS^s[{^)h (B.3) 

Throughout this fit recovery procedure the parameters of V{Nch, Nch) in Eq. ([2]) 
remained fixed. 

Numerically stable results with one-million minimum-bias simulated collisions 
required coarse binning of the Monte Carlo distributions dN(.vt,D / dN^l'^ and 
dNevt,D / d{Npart/'2Y^^ prior to evaluation of Eqs. (IB. ID - (]B.3p . When parameters 
Cpt^m and ainei were shifted excellent fit recovery was achieved with centrality bins 
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Table 5. Same as Table [2] except for Au-Au collisions at ^s^vat = 20 GeV. 







An A 
AVI- A 


n on a\r 








Centrality (%) 




(b) 


( -^part ) 






V 


90-100 


1.9 


1/1 QQ+0-27(l-9) 
I^.OO q 3Q(2.1) 


.+0.1(3.9) 
'-'■^-0.1(1.9) 


9 9+0.1(6.0) 
^■^-0.1(3.3) 


1.29; 


f0.02(1.9) 
-0.02(1.3) 


80-90 


4.8 


iq ^c-+0.27(2.0) 
IO.UJ q 27(2.0) 


^ 0+0-4(6.3) 
^•^-0.4(5.3) 


p. 9+0.4(8.5) 
^^•^-0.3(6.7) 


1.51; 


f0.04(2.4) 
-0.03(1.9) 


70-80 


10.5 


.9 79+0.24(1.9) 
IZ. 27(2.1) 


14 fi+0-9{6.0) 
-^^•"-0.6(4.2) 


.9 ^+1.2(9.1) 
-^'^•°-0.7(5.8) 


1.75; 


f0.06(3.3) 
-0.04(2.1) 


60-70 


20.5 


11 70+0-22(1.9) 
XX. 1 ^_o.26(2.2) 


97 c+l-5(5.5) 
-^'•"-1.1(4.1) 


98 /|+2-6(9.3) 
^°-^-1.9(6.7) 


9 p,^+0.08(4.0) 
^•^"-0.06(3.1) 


50-60 


36.2 


in 09+0-20(1.8) 
iU.OZ_Q 23(2.1) 


.7 0+2.2(4.5) 
^'•'^-1.5(3.2) 


.7 .+4.9(8.5) 
"-^'•^-3.5(6.1) 


2.42; 


f0.10(4.2) 
-0.09(3.5) 


40-50 


59.2 


q 70+0.19(1.9) 

' "-0.20(2.1) 


9+2.6(3.5) 
' "-^"^-1.9(2.5) 




2.83; 


f0.12(4.4) 
-0.11(3.7) 


30-40 


91.6 


1 f\ "1 / "1 f\\ 

^9+0-16(1. 9) 
°-"^-0.18(2.1) 


1 1 q 9+3.1(2.7) 
l±O.Z_2.2(i.9) 


1 Se; 1 +12.3(6.6) 
iOJ.l g 2) 


3.27; 


f0.14(4.2) 
-0.13(3.9) 


20-30 


135.9 


y r,y+0. 15(2.1) 

'■'^'-0.15(2.1) 




qn/l 9+16.7(5.5) 
'JU^-^_14.9(4.9) 


3.72; 


fO. 15(4.0) 
-0.14(3.8) 


10-20 


196.1 


p. ^.+0.11(2.0) 
"-^•"-^-0.12(2.2) 


990 1 +2-8(1.2) 
zzy.x_2. 3(1.0) 


470 5+24.6(5.1) 
^(y.J_22.1(4.6) 


4.19; 


f0.17(4.1) 
-0.16(3.9) 


5-10 


255.4 


Q qR+0-07(1.8) 
"^•^"-0.08(2.1) 


9Q9 9+2.2(0.8) 
^^^•^-2.1(0.7) 


""^•^-30.9(4.7) 


4.55; 


f0.19(4.1) 
-0.18(4.0) 


0-5 


309.2 


9 90+0.07(3.2) 
^•-^^-0.07(3.1) 


q/19 9+2.5(0.7) 
O^Z.Z_2. 7(0.8) 


895 0+38-2(4.6) 
^-38.5(4.7) 


4.83; 


f0.19(4.0) 
-0.20(4.1) 



corresponding to 0-5%, 5-10%, 10-40%, 40-70% and 70-100% total cross section fractions. 
When diffusivity parameter z^i^m was shifted the optimum binning was 0-5%, 5-10%, 10- 
20%, 20-40%, 40-60%, 60-80% and 80-100%. The power-law representation in Eq. (IRTD 
exploits the approximately uniform statistics throughout the domain of the integrands, 
making numerical stability of the results less problematic. 

Variations in the density parameters and Omei produced typical changes in {Nch) 
as a function of centrality of order several percent, up to 5%. Fit recovery in terms 
of (A^c/i) was typically within a few tenths of a percent and always less than about 1% 
where the magnitude of 5{dNch/ drj) relative to dNch/drj was less than about 1%, 4% 
and 1.5% for changes in Cpt^m, Zpt^m and aineu respectively, for all six collision systems. 

A similar method was used to account for shifts in dNevt/dNch representing the 
systematic uncertainties in trigger inefficiency, background contamination, and collision 
vertex reconstruction inefficiency. For these simulations dNi.vt/d{Npart/'^y^'^ remained 
fixed while dN(,^t/dN^[^ varied. The required correction is given by, 

K^)=KKR^s[{^)^-M.,^n^,[{^)^. (B.4) 
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Table 6. Same as Table [2] except for Cu-Cu collisions at y^sjfjj — 200 GeV. 











onn a\r 

Z\J\J Kjrc V 










Centrality (%) 






(b) 


( -^part ) 


{Nbin) 




u 




90-100 


2.2 


10.18 


+0.22(2.1) 

-0.22(2.1) 


^ ^+0.2(7.0) 
-^•"-0.2(7.1) 


. ^.+0.2(10.8) 
-^•'^-0.2(10.6) 


1.18" 


f0.04 
-0.04 


3.1) 
3.0) 


80-90 


4.7 


9.79 


f0.22(2.2) 
-0.20(2.0) 


o 0+0.2(4.5) 
'-'•°-0.2(4.7) 


9 ^+0.2(6.6) 
^•"-0.2(6.7) 


1.35" 


f0.02 
-0.02 


1.8) 
1.8) 


70-80 


8.4 


9.10 


f0.20(2.2) 
-0.19(2.1) 


^+0.2(2.3) 
"■"-0.2(2.8) 


^+0.2(4.1) 
^•"-0.2(4.4) 


1.53" 


f0.03 
-0.03 


1.7) 
1.8) 


60-70 


14.3 


8.40 


f0.19(2.2) 
-0.18(2.1) 


0+0.3(2.6) 
-^"•'^-0.3(2.9) 


Q 0+0.4(4.5) 
^"^-0.5(4.9) 


1.71" 


f0.03 
-0.04 


2.0) 
2.2) 


50-60 


23.3 


7.68 


f0.17(2.2) 
-0.14(1.8) 


. 7 ^+0.4(2.1) 
^ ' -^-0.5(2.6) 


. ^ 7+0.7(4.3) 
-^"•'-0.8(5.1) 


1.94" 


f0.05 
-0.05 


2.4) 
2.5) 


40-50 


36.6 


6.92 


fO. 14(2.0) 
-0.14(2.0) 


.+0.5(1.8) 
'^"■-^-0.6(2.4) 


^°-^-1.5(5.2) 


2.21" 


f0.06 
-0.06 


2.8) 
2.8) 


30-40 


55.4 


6.10 


f 0.1d(2.1) 
-0.13(2.2) 


qo 1 +0.6(1.7) 
-^-0.8(2.0) 


AO 9+2.3(4.8) 
^°-^-2.5(5.1) 


9 p.q+0.08 
^•"-^'^-0.08 


CO OJ 


20-30 


81.5 


5.15 


fO. 10(2.0) 
-0.11(2.2) 


c-q Q+0-8{l-5) 
^"^■^-0.9(1. 7) 


77 q+3.9(5.0) 
' ' -^-3.9(5.0) 


2.89" 


fO.lO 
-0.10 


3.4) 
3.4) 


10-20 


117.4 


3.97 


f0.09(2.2) 
-0.07(1.9) 


7A 4+0-8{l-l) 
'^•^-1.1(1.5) 


190 n+5.6(4.6) 


3.31 


fO.ll 
-0.12 


3.4) 
3.5) 


5-10 


152.2 


2.80 


f0.05(1.9) 
-0.05(1.9) 


QO 9+0.8(0.8) 
^'^•'^-1.0(1.1) 


169 0+7-4(4.3) 
-^"^•^-7.9(4.7) 


f;^+0.13 
'^•"•^-0.13 


3.6) 
3.6) 


0-5 


184.8 


1.75 


f0.11(6.5) 
-0.10(5.9) 


106 5+i-5(^-^) 


9.. .+10.1(4.8) 
-^J-J--+_10.5(5.0) 


3.97 


f0.15 
-0.15 


3.8) 
3.8) 



Appendix C. 

The Monte Carlo simulation model discussed in Sec. [H] for transverse particle production 
and energy deposition in the STAR CTB for p-p and Au-Au minimum bias collisions 
is described here. For p-p collisions the interaction point along the beam-line was 
randomly selected within ±25 cm of the geometrical center of the STAR detector |44j . 
The minimum-bias charged particle multiplicity was obtained by sampling the negative 
binomial distribution (NBD) fitted to the UA5 = 200 GeV p-p data ^ [H] 
with parameters (n) = 21.6 and k = 4.6. The pt, r] and particle species (tt^, 
K^, proton or antiproton) for each charged particle produced in the collision was 
sampled from the measured dNch/dpt, dNch/drj, and particle species distributions, 
respectively. Distribution dNch/dpt was assumed to be proportional to exp(— /?mt), 
where mt = \/ pi + m^, m is the particle rest mass, and inverse effective temperature 
fi = 5.7, 5.2 and 4.8 GeV~^, for pions, kaons and protons, respectively. The P values 
were obtained by fitting data [IS] within the approximate range irit — m < 0.6 GeV. 
Pt was restricted to be less than 2 GeV/c. Distribution dNch/drj was obtained from 
the UA5 [18] measurements for y/s = 200 GeV non-singly diffractive p-p collisions 
for |?7| < 4.6. The particle species probabilities for pions, kaons and protons at low 
momentum were assumed to be 0.85, 0.085 and 0.065, respectively [23l B6] . 

Particle trajectories were extended outward from the collision point and 



Applicability of Monte Carlo Glauber models to relativistic heavy ion collision data 23 

Table 7. Same as Table [2] except for Cu-Cu collisions at y^sjfjj = 62 GeV. 
Cu-Cu 62 GeV 



Centrality (%) 




(b) 




{Nbin) 




V 


90-100 


1.8 


1 no ^ ^ 

IU.UU q 21(2.1) 


-^0 1 ("^ ZL^ 

^•°-0.1(2.7) 


' -0.1(4.4) 


1 


4-0 09('1 ZL^ 

^^-0.02(1.3) 


80-90 


3.8 


^•"^-0.20(2.1) 


'^•°-0.1(3.3) 


„ „4-0 1 (d 2) 
^•"-0.1(4.1) 


1 


„ , 4-0 02f 1 1 "I 
-0.02(1.3) 


70-80 


6.7 


„ „ —-1-0 21 (2 S"! 
"•"^ ' -0.17(1.9) 


„ -4-0 1 (2 2) 
"•^-0.2(3.2) 


. „4-0 2f3 3~) 
^•^-0.2(4.7) 


1 


_ „4_n 02f 1 
'-'^-0.02(1.6) 


60-70 


11.2 


„ -4-0 1 7(2 1 "1 
-0.16(2.0) 


, „ „4-n 3('2 7'! 
-^"•"-0.3(2.8) 


„ „4-0 4^4 H) 
O-O_0.4(4.5) 


1 


""-0.03(1.8) 


50-60 


17.9 


4-0 1 

' -0.15(2.0) 


, Lf) d(2 9~t 

1 c y-l-u.'uz.zj 

' -0.5(2.9) 


-t-D 7(A 

-^^•^-0.8(5.2) 


1 


of-+0. 04(2.3) 
"^-0.04(2.4) 


40-50 


27.7 


"•°^-0.13(1.9) 


oc 0+0.4(1.8) 
^^•'-'-0.7(2.6) 


9fi q+l-l{4.3) 
^"•'-'-1.4(5.2) 


2 


p,Q+0.05(2.6) 
^"-0.06(2.7) 


30-40 


41.3 


"•^^-0.11(1.8) 


q^ ^ + 0.5(1.5) 
'^"•°-0.8(2.2) 


^-^1.9(4.3) 
^'-'•^-2.2(5.1) 


2 


q p. -1-0. 07(3.0) 
"^-0.07(3.1) 


20-30 


60.0 


p. „p+0.11(2.1) 
^•^"-0.08(1.6) 


P.^ .-h0.5(1.0) 
"^-^•-^-1.0(1.9) 




2 


1^ p. -f 0.07(2. 7) 
""^-0. 09(3. 3) 


10-20 


84.9 


o +0.07(1.8) 
"^•^-^-0.06(1.6) 


7. ^-h0.7(1.0) 
' -^-"-l. 0(1.4) 




2 


qy+0. 09(3.1) 
^ -0.11(3.6) 


5-10 


108.9 


^ 7^-f0.05(1.7) 
'"-0.03(1.3) 


on i+0-8{0-9) 

^•J- -^-1.1(1.2) 


14c; 7+6.5(4.5) 
-^^'-'■'-6.9(4.7) 


q r)q + 0.12(3.6) 

"^•^"^-0.12(3.6) 


0-5 


132.2 


„„+0. 10(5.6) 
^- ' ' -0.08(4.7) 


inQ 9+l-4(l-3) 
-^'^'^•'^-1.5(1.5) 


17^-4-9.l(5.1) 


q .^-1-0.14(4.0) 
'^•^'-'-0.14(3.9) 



Table 8. Absolute values ol individual error contributions (in percent) to collision 
geometry bin average quantities for Au-Au collisions at ^/snn — 200 GeV relative to 
the reference values in Table [5] as explained in the text. Errors are denoted as A (6), 
etc. The left-most column lists the sources of uncertainty. Average errors for centrality 
bins 60-100%, 20-60% and 0-20% are listed from left to right, respectively, for each 
instance. 



Error 




Au-Au 200 GeV Errors 




Source 


A(6)(%) 


^{Npart){%) 


A(iVM„)(%) 


Az.(%) 




0.9, 1.0, 0.8 


0.1, 0.1, 0.2 


0.1, 0.9, 2.2 


0.0, 0.9, 2.0 




1.7, 1.6, 1.5 


3.2, 1.6, 0.3 


5.2, 5.1, 2.6 


2.1, 3.5, 2.3 




0.1, 0.1, 0.0 


0.4, 0.5, 0.3 


1.0, 2.0, 2.5 


0.6, 1.5, 2.2 


X 


0.0, 0.0, 0.1 


0.2, 0.0, 0.0 


0.3, 0.0, 0.0 


0.2, 0.1, 0.0 


a 


0.0, 0.0, 0.2 


0.3, 0.1, 0.0 


0.4, 0.2, 0.1 


0.1, 0.1, 0.0 


V-iovva 


0.1, 0.0, 0.7 


2.2, 0.1, 0.2 


3.5, 0.3, 0.3 


1.2, 0.2, 0.1 


trigger 


0.0, 0.0, 0.0 


0.1, 0.0, 0.0 


0.1, 0.0, 0.0 


0.0, 0.0, 0.0 


vertex 


0.0, 0.0, 0.1 


0.1, 0.0, 0.0 


0.1, 0.0, 0.0 


0.0, 0.0, 0.0 


tracking 


0.0, 0.0, 0.1 


0.0, 0.0, 0.0 


0.0, 0.0, 0.0 


0.0, 0.0, 0.0 



approximately transverse to the beam direction toward the STAR Central Trigger Barrel 
(CTB) detector [391 Hlj which was approximated by a uniform cylinder 220 cm in radius, 
coaxial with the beam line, and having full 2ii azimuthal coverage and total longitudinal 
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Table 9. Same as Table |8] except absolute values of errors (in percent) for Cu-Cu 
collisions at y^sJ^J^ — 200 GeV relative to the reference values in Table [H] as explained 
in the text. 



Error 




Cu-Cu 200 GeV Errors 




Source 


A(&)(%) 


A(A^part)(%) 


A(7Vfc„)(%) 






1.0, 1.0, 1.0 


0.2, 0.4, 0.5 


0.2, 1.1, 2.6 


0.2, 0.8, 2.0 




1.9, 1.8, 1.7 


2.1, 1.8, 0.8 


3.2, 4.3, 3.2 


1.2, 2.5, 2.3 




0.2, 0.2, 0.2 


0.2, 0.4, 0.4 


0.5, 1.5, 2.2 


0.3, 1.1, 1.8 


X 


0.0, 0.0, 0.0 


0.0, 0.0, 0.0 


0.2, 0.0, 0.1 


0.2, 0.1, 0.1 


a 


0.1, 0.1, 0.3 


0.8, 0.2, 0.1 


1.2, 0.4, 0.2 


0.4, 0.2, 0.1 


7^-form 


0.3, 0.1, 1.9 


2.8, 0.1, 0.5 


4.7, 0.5, 0.6 


1.5, 0.5, 0.3 


trigger 


0.0, 0.0, 0.1 


0.2, 0.1, 0.0 


0.2, 0.1, 0.0 


0.1, 0.0, 0.0 


vertex 


0.0, 0.0, 0.0 


0.0, 0.0, 0.0 


0.2, 0.0, 0.0 


0.0, 0.0, 0.0 


tracking 


0.0, 0.0, 0.0 


0.0, 0.0, 0.0 


0.0, 0.0, 0.0 


0.0, 0.0, 0.0 



Table 10. Percent changes in collision geometry bin average quantities relative to the 
reference values in Tables [2] - [7] due to background contamination of the peripheral 
collision events as explained in Sec. SI Background contamination errors for centrality 
bins not listed were negligible. 



Background Contamination Errors 

System Centrality A(&> l::^{Npart) ^{Nun) Az/ 





(percent) 


(%) 


(%) 


(%) 


(%) 


Au-Au 200 GeV 


90-100 


-0.3 


5.6 


7.9 


1.8 


23% contamination 


80-90 


0.2 


-1.6 


-1.7 


-0.1 


in 88-100% bin 












Au-Au 130 GeV 


90-100 


-0.3 


4.8 


6.7 


1.6 


22% contamination 


80-90 


0.2 


-1.5 


-1.5 


0.0 


in 91-100% bin 












Au-Au 62 GeV 


90-100 


-0.1 


1.3 


1.6 


0.3 


18% contamination 


10-50^ 


-0.1 


0.2 


0.4 


0.1 


in 96-100% bin 












Au-Au 20 GeV*^ 


90-100 


-0.1 


2.7 


3.5 


0.8 




80-90 


0.0 


0.0 


0.1 


0.0 




40-70^ 


-0.1 


0.5 


0.7 


0.2 


Cu-Cu 200 GeV*^ 


90-100 


-0.2 


1.7 


2.3 


0.8 


Cu-Cu 62 GeV'^ 


90-100 


-0.1 


1.5 


2.0 


0.6 



Average percent changes within the combined centrality bin. 



^ Reference uncertainties assuming 10% background contamination in the number of 
collisions in the nominal 90-100% centrality bin, equivalent to 1% overall background. 



length of 484 cm located symmetrically about the geometrical center of the STAR Time 
Projection Chamber [33]. Charged particles were propagated in a 0.5 T solenoidal 
magnetic field as helices until they either intersected or missed the CTB. The simulated 
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Table 11. Upper half-max end-point positions Npart,max and Nun. max from the 
power-law distributions da/d{Npart/'2)^^'^ and da/dN^I^, respectively, for the six 



collision systems 


studied here. 




System 


A^ part,max 


^ bin,max 


Au-Au 200 GeV 


379.3 


1166 


Au-Au 130 GeV 


377.5 


1082 


Au-Au 62 GeV 


374.9 


988 


Au-Au 20 GeV 


373.4 


926 


Cu-Cu 200 GeV 


115.8 


230 


Cu-Cu 62 GeV 


112.7 


196 



mean number of charged particles per unit pseudorapidity at rj = intersecting the 
CTB was consistent with the measured Upp from UA5 [18j. 

Energy deposition in the plastic scintillator material (CH) was estimated by 
sampling the Landau distribution [37] for the calculated most-probable energy loss [H] 
for 1 cm thick plastic taking into account the range of crossing angles between the helical 
trajectories and the cylindrical detector. The areal density (effective thickness) along the 
path through the material was estimated to be (1.032 gm/cm^)/ cos(a) [16] where a is 
the angle of incidence between the particle trajectory helix and a normal to the detector 
cylinder at the intersection point. The mean electron excitation (ionization) energies 
for carbon and hydrogen were assumed to be 80 eV and 20 eV, respectively [I6|. Energy 
depositions in the carbon and hydrogen components of the detector material were added 
together. Light attentuation in the plastic scintillator (380 cm attenuation length [39] ) 
was included in the simulated CTB output signal. For each p-p colhsion the total output 
from each produced charged particle which deposited energy in the CTB scintillators 
was summed, resulting in a simulated result for the quantity SUM in Sec. |6l The p-p 
results in Fig. H] correspond to one-million triggered events with non-zero CTB signal. 
The resulting trigger output frequency distribution, dNtrig{pp) / dSUM , served as input 
for the Au-Au trigger simulations described next. 

For Au-Au collisions the Monte Carlo Glauber model presented here was used to 
generate an ensemble of ^/SNN = 200 GeV simulated minimum-bias collisions. For 
each Au-Au collision the dNtrigipp) / dSUM distribution was sampled Nch/npp times 
corresponding to the average number of p-p collisions required to generate mean 
multiplicity N^h- The sampled values of SUM were added to obtain the total CTB 
trigger output for each simulated Au-Au collision. The results in Fig. |4] correspond to 
one-million Au-Au triggers with non-zero CTB output. 
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